Thermal Conductivity of Isotopically Enriched Si: Revisited 
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The thermal conductivity of isotopically enriched Si (enrichment better than 99.9%) 
was redetermined independently in three laboratories by high precision experiments on 
a total of 4 samples of different shape and degree of isotope enrichment in the range 
from 5 to 300 K with particular emphasis on the range near room temperature. The 

results obtained in the different laboratories are in good agreement with each other. 
They indicate that at room temperature the thermal conductivity of isotopically 
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enriched Si exceeds the thermal conductivity of Si with a natural, unmodified isotope 
mixture by 10±2 %. This finding is in disagreement with an earlier report by Ruf et al. 
At -26 K the thermal conductivity of Si reaches a maximum. The maximum value 
depends on sample shape and the degree of isotope enrichment and exceeds the thermal 
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conductivity of natural Si by a factor of ~8 for a 99.982% Si enriched sample. The 
thermal conductivity of Si with natural isotope composition is consistently found to be 
-3% lower than the values recommended in the literature. 
PACS: 66.70.+f, 63.20.Mt, 74.25.Fy, 65.40.Ba, 65.90.+i 
Keywords: silicon, stable isotopes, thermal conductivity 

INTRODUCTION 

Phonon scattering due to the presence of different isotopes in an otherwise pure crystal (no 
chemical defects or dislocations) has been identified as a mechanism that strongly affects the 
thermal conductivity K [1-3]. The availability of larger quantities of highly isotope-enriched 
materials recently revived the interest in this effect in order to study the mechanisms underlying 
the thermal conductivity and possibly to improve material properties [4-8]. 
Recently, the observation of a significantly enhanced thermal conductivity of isotopically 
enriched Si near room temperature by Capinski [9] and Ruf et al. [10] has generated large 
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interest. In these studies, it was found that the thermal conductivity of isotopically enriched Si 
is enhanced by about 60% over that of silicon with natural isotopic composition. This rather 
high isotope effect attracted considerable attention concerning both, fundamental physics and 
applications. For technical applications, a significantly enlarged thermal conductivity at room 
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temperature would be of interest for high performance electronic devices. [11]. From the 
fundamental physics aspect, the experimental results were unexpected since the observed 
isotopic effect was significantly larger than the prediction of simple theoretical estimates 
[10-13] and more advanced model calculations [14]. On the other hand, theoretical papers 
[15-17] were published with the results supporting the data of refs. [9,10]. A large isotopic 
effect at room temperature is in principle only possible if normal phonon-phonon scattering 
processes play an important role in determining the formation of the non-equilibrium 
distribution function of phonons at these temperatures. 

The large room temperature isotope effect was questioned following the results of an 
experimental study by Gusev et al. [18] which indicate that the room temperature isotopic effect 
amounts to only 7%, and more recently by an erratum in which an increase at room temperature 
of about 10% was reported [19]. Thus, presently there is considerable uncertainty concerning the 
effect of isotope disorder on the thermal conductivity of silicon, especially in the technologically 
relevant regime around room temperature. This and the broad interest for potential technical 
applications led us to carefully redetermine the thermal conductivity of Si in comparison with 
that of Si with the natural isotope composition ("^'Si). For this purpose, absolute measurements 
of the temperature dependence of the thermal conductivity of several samples of Si and of 
natural silicon using different experimental techniques, were carried out independently in three 
laboratories, in Stuttgart, Nizhny Novgorod and Moscow. Here we report the results of these 
experiments. All experiments conclusively find that at room temperature the thermal 
conductivity of ^^Si exceeds that of "^'Si by only 10% while at the maximum close to 26 K, icis 
enhanced by almost an order of magnitude. 

EXPERIMENTAL DETAILS 

Thermal conductivities were measured on a total of four different bar-shaped samples of 
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enriched Si, (99.9% and higher) with a broad range of cross sections (-4-20 mm ) and 
lengths (20 - 50 mm). The samples were cut from crystals grown by the floating-zone method 
such that the orientation along the sample length was always [100]. A detailed description of the 
preparation of isotopically enriched Si crystals has been given elsewhere [20]. The samples 
with natural isotopic composition, which were measured for comparison, had the same sizes and 
crystallographic orientation and similar chemical purity. Chemical, isotope and geometrical 
details of all investigated samples of ^*^Si are compiled in Table I. The surfaces of samples ^*^Si- 
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M and ^'Si-NN were particularly ground with 14 |im abrasive powder slurry to ensure diffuse 
scattering of thermal phonons from the sample surfaces at low temperatures. 
In all three laboratories, the thermal conductivities were measured by the steady-state heat flow 
technique using either platinum resistance thermometers (Lake Shore Cryotronics, Inc.) (50- 
320 K) or Cemox thermometers (CX-1050 Lake Shore) (5-300 K) as temperature sensors and a 
surface mounted thick film (SMD) resistor as a heater which was glued to one end face of the 
samples. The opposite end of the sample bars was thermally anchored to a Cu block which was 
thermally coupled to the cryostat. Distances between the temperature sensors were Ax=10-20 
mm. After sufficient stabilization to a particular temperature T, the electrical power P dissipated 
in the SMD heater was adjusted such that typical temperature gradients Ar amounted to 0.01- 
0.2 K (Kurchatov Institute and IChHPS RAS) while in Stuttgart several (typically 4-6) stable 
gradients 0.01 < AT < 0.7 K were generated and the thermal conductance was obtained from a 
least-squares fit of the set of measured {Pi, ArO pairs. The latter approach provides additional 
signal-to-noise improvement and enables a critical assessment of deviations from the linear 
relationship ATi cx Pj and thus allows us to recognize self -heating of the sample due to possible 
insufficient thermal coupling to the bath. The thermal conductivity KiT) at an average 
temperature (Ti+T2)/2 in turn is calculated with the equation P = k(T) A AT/ Ax, where AT = (Ti 
- T2) is the stable temperature gradient measured between the two thermometers at a distance Ax 
apart in the presence of the constant heat flow P= U -I, I being the current and U the voltage drop 
across the SMD resistance heater. A is the cross section of the sample. The samples were 
mounted into evacuated Cu cans which were immersed into either liquid nitrogen or liquid He 
(IChHPS RAS, Kurchatov Institute) or mounted into a continuous flow He evaporation cryostat 
(MPI FKF). 

The temperature sensors were either attached to the sample bars by Cu clamps with an indium 
layer to improve the thermal contact (Kurchatov Institute and IChHPS RAS) or by glueing 
either the Cemox or Pt chip thermometers directly to the sample using an epoxy resin (MPI 
FKF). The absolute errors in the determination of the thermal conductivity arise from the finite 
thickness of the temperature sensors or the clamps employed to connect the sensors to the 
samples. For both methods we estimate the absolute errors to 2-3%. Particular emphasis was 
given to using the same thermometer mounting when measuring the isotopically enriched and 
the natural Si samples. This enables us to reduce the relative error by a factor of -2 when 
comparing both types of samples. Towards higher temperatures, thermal radiation opens an 
additional path to dissipate the supplied heat. Thermal radiation increases with temperature like 
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T ^Ar and, compared to thermal conduction, it becomes increasingly important for samples with 
smaller cross sections and comparatively larger surface areas. For the sample with the smallest 
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cross section investigated (MPI FKF, Si-Sl, A~3.8 mm ) assuming an infrared emissivity of 
-0.5, radiation losses were estimated to contribute by about 2% to the heat dissipation at room 
temperature. Thermal radiation losses at room temperature are therefore smaller than other 
experimental errors and become negligible for the samples with larger cross sections. To test the 
resolution and eventually correct for radiation losses, the thermal conductivity of isotopically 
enriched samples with different cross sections was determined with especially high precision 
around room temperature and above. 
RESULTS AND DISCUSSION 

The thermal conductivities k(T) of three isotopically enriched ^^Si samples as measured 
independently by the three laboratories are displayed in Fig. 1 in comparison to the thermal 
conductivities of Si with the natural isotope composition and to the literature data [21]. For all 
sets of data, the thermal conductivity of ^^Si at 300 K is 10±2% larger than that of °^'Si. This 
finding is in very good agreement with the results of Morelli et al. who predicted an increase at 
room temperature of 12% [14]. Above 80 K, the results of the different laboratories agree to 
within 3%. At lower temperatures, we observe characteristic differences, which can be ascribed 
to the different degrees of isotope enrichment, chemical purity and, at the lowest temperatures, 
to the different sample dimensions and surface finish. Comparison of the thermal conductivity of 
"^*Si with the standard data recommended by Touloukian [21] additionally reveals that at room 
temperature our three laboratories consistently find the thermal conductivity of "^'Si to be about 
3% smaller than the standard values. Figure 2 displays the inverse thermal conductivity 
('thermal resistance') of "^*Si together with a polynomial of order two as guide for the eyes. In 
the temperature range 220 K < T <320 K the thermal resistance, within error bars, closely 
follows the equation K '\m K W"^] = l.60(2)x10"^ T + 2.37(6)x10"^ T ^ with coefficients 
similar to those derived by Glasbrenner and Slack [17], however, truncating a temperature 
independent contribution of -3x10"^ mK W"^ as used in [13]. 

The low-temperature data displayed in Fig. 3 reflect the degree of isotope enrichment, chemical 
purity and surface conditions of the different samples under investigation. At low temperatures, 
the thermal conductivity approaches a T power law behavior, however, with slightly different 
prefactors for the different samples. While at these temperatures the thermal conductivity of the 
^^Si sample investigated at the Kurchatov Institute and at the IChHPS RAS coincides with that 
of "^*Si, the thermal conductivities of the MPI FKF ^^Si sample are by about a factor of 2 larger 
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than those of ""'Si. We ascribe this difference to additional phonon scattering by chemical 
impurities or different surface scattering due to different sample surface finishing [7]. We finally 
discuss the thermal radiation losses which around room temperature may provide an additional 
channel for heat dissipation in samples with a small cross section. Figure 4 displays the 
"nominal" thermal conductivities of two samples which differ in cross section by about a factor 
of 5. The total sample length and the distance between the temperature sensors were 
approximately the same. We can clearly observe an increasing difference of the "nominal" 
thermal conductivity between the two samples which grows approximately ocT^ thus indicating 
that the difference is due to increased radiation losses for the sample with smaller cross section. 
Using the Stefan-Boltzmann law for thermal radiation, we can estimate the ratio of heat flow due 
to thermal radiation Prad and conduction Pcond to be P^^^ / P^^^^ = a £ IS /( A k(T )), where S 

is the sample surface, A the cross section, / the sample length, ^ the infrared emissivity and a the 
Stefan-Boltzmann radiation constant. Using the sample dimensions and assuming an average 
infrared emissivity e~Q.5, we estimate that for the sample with cross section A~4 mm the 
radiation losses at room temperature amount to about 2%, which is in very good agreement with 
the experimental observation. Consequently, for the larger sample of (A~20 mm ) radiation 
losses at room temperature are negligible. 

In summary, in a collaboration involving three laboratories, we have independently 
redetermined the thermal conductivity of natural Si and isotopically enriched Si samples. In 
the temperature range 80 K<r<300 K, we find excellent agreement between the three sets of 
results. We therefore conclude that at room temperature the thermal conductivity K of Si 
exceeds Koi "^'Si by only 10±2 %. Optimized experimental conditions, improved samples and 
the close coincidence of all experimental results indicate that the original measurements by Ruf 
et al. represent an overestimate of k(T) at room temperature, a fact which has already been 
admitted by these authors [19]. In addition, we find that close to room temperature the values of 
the thermal conductivity of "^'Si are about 3% lower than those recommended in the existing 
literature. 
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Fig. 1 Thermal conductivities (o) of three samples of Si. (a) sample Si-M, (b) sample 
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Si-NN and (c) sample Si-S2 (sample details cf. Table 1) in comparison with the thermal 
conductivities of samples of "'"Si with similar sample sizes (~). The solid lines represent the 
thermal conductivity of "'"Si given in [21]. 
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Fig. 2 Themial resistivities of three samples of Si ( o Si-S2, • Si-NN, p Si-M) in comparison with 
the thermal resistivity of "^'Si (■,□). The full lines are fits with 2"'^ order polynomials in the temperature 
range 220 K<r<320 K with parameters: ^^Si: ic'^ [m W K"^]= 1.33(2)xl0"^ T + 2.66(5)xl0"^ T \ "''Si: 
K'^ [mW K"^]=1.60(2)xl0"^ 7+ 2.37(6)xl0"^ T^. The dashed line represents the standard data for ""'Si 
compiled in [21]. 
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Fig. 3 Thermal conductivity of two samples of Si (o Si-M, • Si-S2) in comparison with 
the thermal conductivity of "^'Si with size and surface finish similar to sample ^^Si-M. 
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Fig. 4 Thermal conductivity of two samples of ^^Si with cross section 3.8 mm^ (A, sample ^^Si-Sl) and 
19.6 mm^ (o, sample ^^Si-S2). 
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isotopic composition (%) 




content of impurities 


temperature range 
(K) 


size 


sample 


^^Si 


^'^Si 


^•^Si 


orien- 
tation 


carbon, 
oxygen 
(cm ) 


metals 
(cm ) 


(mm^) 


^^Si-M 


99.9829 


0.0144 


0.0027 


[100] 


2x10^^ -C 
1x10^^ -O 


<10"^ 


5 - 300 


2.5 X 2.0 x 23 


^^Si-Sl 


99.859 


-0.13 


-0.02 


[100] 


5x10^^ -C 
2x10^^-0 


<10"^ 


50 - 320 


2.0 X 1.9x60 


^^Si-S2 


99.979 


-0.019 


-0.002 


[100] 


5x10'^' -C 
2x10^^-0 


<10"^ 


6-320 


4.5 X 4.35 X 60 


^^Si-NN 


99.9829 


0.0144 


0.0027 


[100] 


2x10^^ -C 
1x10^^ -O 


<10"^ 


5-300 


3.12x2.0x20 



Table 1 Physical and chemical parameters of the Si samples measured at RRC, Kurchatov Institute (Moscow) ( 
^^Si-M), at the MPI FKF (Stuttgart)( ^^Si-Sl, ^^Si-S2) and at the IChHPS RAS (Nizhny Novgorod) ( ^^Si-NN). 



